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Abstract 

The systems with multimode nonstationary Hamiltonians quadratic in position and momentum 
operators are reviewed. The tomographic probability distributions (tomograms) for the Fock states 
O \ and Gaussian states of the quadratic systems are discussed. The tomograms for the Fock states are 
expressed in terms of multivariable Hermite polynomials. Using the obvious physical relations some 
O ■ new formulas for multivariable Hermite polynomials are found. Examples of oscillator and charge 
£^ \ moving in electromagnètic field are presented. 

VO : 

°i Introduction 

^The tomographic probability distribution was introduced Q , for reconstructing the Wigner function || 
çrof quantum states. The optical tomography scheme Q was used in experiments for measuring quantum 
states of photons. Recently the symplectic tomography method was suggested ||. The possibility to 
^jjlescribe the quantum state in terms of tomographic probability distribution was employed to reformulate 
^tRiantum mechanics f| , avoiding the complex probability amplitudes, wave function and density matrix 
conceptual ingredients. The evolution equation for tomograms ||, |7J and equation for energy levels 
*^*n! terms of the tomograms |§ is tomographic counterpart of Moyal equations |§ written for the Wigner 
^f unction of quantum system. The tomographic probability distribution was introduced also for systems 
Pyith spin degrees of freedom ( []I0] , [pH , |Ï2| , |Ï3| , [[14]]). The discrete spin degrees of freedom and the 



íontinuous ones like a position were considered in tomographic representation in |Ï5|. The tomograms are 



Güfeed for description of quantum states because they contain the same information on the state as other 
^eharacteristics like Wigner function ||, Glauber-Sudarshan P-function |Ï6| and | Ï7| , Husimi Q-function 



but the tomograms have specific property of Standard probability distribution. Important role play 
?4he quantum systems with Hamiltonians which are quadratic in position and momentum |Ï9| . The charge 



moving in homogeneous elèctric and magnètic field and linear vibrations of polyatomic molècules are 
examples of such systems. The coherent states and Fock states for such systems were studied using time- 
dependent integrals of motion linear in position and momentum in |ÏÏJ. The wave functions of coherent 
states and Fock states were obtained for the quadratic systems in explicit form. On the other hand the 
tomographic description of the states of the multimode quadratic system has not been presented till now. 

The aim of our work is to obtain the tomograms of the specific quantum states of the quadratic 
systems. We find the explicit Gaussian tomograms for coherent (and squeezed) multimode vibrations. 
Also we show that the tomograms of Fock states are the positive probability distributions expressed in 
terms of multivariable Hermite polynomials. The expression for transition probabilities obvious from 
physical point of view generate some new formula for the Hermite polynomials. 

The paper is organized as follows. In Section 2 we present a review of properties of coherent and Fock 
states for multimode systems with nonstationary quadratic Hamiltonians. In Section 3 we discuss the 
coherent and Fock states for such systems in the tomographic probability representation. In Section 4 we 
obtain some new formulas for multivariable Hermite polynomials. The general results are illustrated by two 
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examples in Section 5 where a driven parametric oscillator and a charged particle moving in nonstationary 
elèctric field are considered in tomographic probability representation. Concluding remarks on properties 
of tomographic probability representation for quadratic quantum systems are presented in Section 6. The 
mathematical properties of symplectic matrices are described in Appendix A and calculation of the matrix 
exponent which is used in main text is presented in Appendix B. 



2 Coherent and Fock states of systems with quadratic Hamil- 
tonians 

In this section we review the problem of finding linear integrals of motion and constructing the coherent 
and Fock states of nonstationary iV-dimensional systems described by Schròdinger equation: 



ift^^ = J ff(í)V(*,í), 

with an arbitrary time- dependent quadratic Hermitean Hamiltonian, i.e. 



(1) 



- - 2N 2N 

H(t) = -q T B(t)q + c T (t)q = - q a B aP {t)% + ^c 7 (í)ç 7 



(2) 



a,P=l 



7=1 



where we use 2 N- vector q = (p,x), N- vector p = ,Pn), N- vector x = (xi, , ïjv), and p m = 

—ihd/dx m , x m = x m in position representation, m — 1, N. The notation q is used for column vector, 
and the notation q T is used for row vector which is transposed column vector. The matrix elements B a p(t) 
form the 4-block matrix 



B(t) 



Bpp(t) B px (t) 
B xp (t) B xx (t) 



(3) 



which is a real 2N x 2iV-symmetric matrix, and 2iV-vector c(t) = (c p (t) , c x (t)) with iV-vectors c p (t) = 

(cpi(í), , c P N{t)), c x (t) = (c x i(í), , c X N(t)) are arbitrary real vectors. Coherent and Fock states of the 

system with Hamiltonian (§) can be obtained in the framework of method of time-dependent invariants 
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Let us consider the following operators: 

Pm(t) =Ü(t)p m Ü- 1 (t), 

X m (t) = Ü(t)x m Ü- 1 (t), 



Q a (t) = Ü(t)q a Ü- 1 (t), 

Q(t) = (Qi(t), ,&*(*)) = (P(t),X(t)). 



(4) 
(5) 



where U{t) is unitary evolution operator, which connects wave function ip(x, t) given at the moment t with 
the wave function ip{x, 0) given at the moment t = 



ip(x,t) = U(i)i)(x,0). 
One can verify that the operators Q a {t) satisfy the following equation 

dQ a (t) 



(6) 



ih- 



dt 



H(t),Q a (t) 



(7) 
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It means that the operators Q a (t) are quantum integrals of motion, i.e. their mean vàlues (Q a (t) 
remain constant on an arbitrary solution if) of Schòdinger equation 

d -^Qa(t)) =0. ( 



at \ ' ' / v 

One can verify that Q(t) satisfying (0), where H(t) is a quadratic Hamiltonian, is given by the formula 

Q(t)=A(t)q + A(t), (9) 
where real 2N x 2A-matrix A(í) and real 2N- vector A(t) read 

A pp (í) A px (t) 



A(t) 



Axp(í) A ra (í) 



A(í) 



(10) 



Indeed, substituting this expression into (0), we have 



ihÀ(t)q + ihÀ(t) 



lq T B(t)q + c T (t)q,A(t)q + A(t) 



-A(t) \q T B(t)q,q] + A(t) [c T (t)q,q] , (11) 



where we introduce the following notations ( a and b are arbitrary N- vectors): 



a, b 



\a N J 



a 1 ,b 1 



a 1; 6jv 
ün, bN 



i.e 



a, 6 



(12) 



Using commutation relations for operators q: 






—En 


E N 






(13) 



we can calculate commutators on the right hand side of equation ([□])• Finally we obtain the equation 

À(t)q + A(í) = A(t)E 2N B(t)q + A(t)E 2N c(t). (14) 

From this equation we can conclude that Q(t), given by the expression (§), satisfies (0) provided the 
matrix A(í) and the vector A(t) satisfy the following evolution equations 



A(t)=A(t)E m B(t), 
A(í) = A{t)E 2N c{t). 

The initial conditions for these evolution equations are taken in the form 

A(0) = E 2N , A(0) = 0. 



(15) 



(16) 



The above initial conditions correspond to initial vàlues of the integrals of motion Q(0) = q. The coherent 
states of the quadratic multimode systems play the special role. These states have the Gaussian form and 
they are close to classical states of the vibrating oscillators. The states are labeled by continuous complex 
quantum numbers. The coherent states form the overcomplete nonorthogonal basis in the Hilbert space. 
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To find the coherent states of the system with the Hamiltonian (H), let us introduce the annihilation and 
creation operators 

a = A p p + A x x, à + = A*p + A* x x, (17) 

where A p and A x are time-independent A" x AT-matrix, and 2N annihilation and creation operators a and 
a + satisfy the following commutation relations: 



a \ la 

a+J'U + 



J 2N • 



To satisfy the condition (|ïq) , the matrices A p and A x must possess the following properties: 

A x A T p - A p A T x = 0, (19) 
A X A; - A p At = ~E N . (20) 

In Appendix A we show that such matrices also possess the following properties 

(i) Matrices A p and A x are non-singular, 

(ii) A T V A* = A+Ap, A T X A* X = A+A x , 

(iii) A x A p — A X A* = ApA x — ApA x = ^E N . 

Below we use this properties to transform some expressions. Making transformation — > Ak(t) = 
U(t)aiJJ(t)~ l on both sides of (^) we obtain integrals of motion in explicit form 

À{t) = AQ{t) = Q{t)q + S{t), Q{t) = AA{t), 6{t) = AA{t), (21) 

where the rectangular matrix 

A=\\A P ,A X \\ (22) 
contains time-independent A^ x A^-blocks and the A^-vector 

A{t) = (A 1 (t),...,À N (tj) (23) 

takes the initial value 

2(0) = a. (24) 

Making the same procedure as we make obtaining the equations for the matrix A(í) we obtain the equations 
for the matrix Q(t) and vector S(t) that are similar to (1Ï5|) 

d(t) = n(t)E 2N B(t), 

S{t) = Ü{t)E 2N c{t). 

But the initial conditions for the matrix Q(t) and the vector S(t) are different from (|T6|), i.e. 

íí(0) = A, 5(0) = 0. (26) 
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Let us introduce the following matrices: 

A p {t) = A p A pp (t) + A x A xp (t), 
A x (t) = A p A px (t) + A x A xx (t). 

The matrices A p (t) and A x (t) possess the properties similar to corresponding ones for the matrices A p 
and A x (see Appendix A). Using the matrices we can express the integrals of motion A(t) through the 
momentum and position operators p and x: 

À(t)=A p (t)p + A x (t)x + S(t), (28) 
n(t) = \\A p (t),A x (t)\\. (29) 

Let us introduce the definition of coherent and Fock states of the system with Hamiltonian H(t). For 
all a. = («!, ... , ajv) G there exists the normalised eigenfunction i[) OL (x, t) of operators A k (t) 

A k (t)ip a (x, t) = a k íp a (x, t), (30) 

which is called the wave function of the coherent state; and for all n = (jii, ... , n/v) G N w there exists 
the normalised eigenfunction ip n (x,t) of operators A^{t)A k {t) 

À+(t)À k (t)ip n (x, t) = n k ^j n (x, t), (31) 

which is called the wave function of the Fock state. It is well known that we can represent the coherent 
state in the form of series: 

°° ct m 

ip a (x,t) = e-^ a ll y2if, m (x,t)-=. (32) 

' ' a / i-n 



m=0 



Here we use the following notations: 

N N N 



33! 

i=l i=l i=l 



It is easy to show that for all a = (ai, ... , a/v) G C N the function f a (x, t), defined as 

Ux,t) = ip ct (x,t)e^\ (34) 
is an analytic function of complex argument a. and it satisfies the following equations: 

À k (t)Ux,t)=a k f a (x,t), À+(t)f a (x,t) = dfa ^ f \ (35) 

oa k 

Indeed, we can represent the function f a (x,t) in the form of the series 

2 °° ct m 
f a (x,t) = ifj a (x,t)eï M = Y] i) m (x,t)—=. (36) 

n ym\ 

m=0 v 

It follows from this expression that the function f a (x,t) is analytic function. To prové the validity of 
the action of the operators A m (t) and on the function f a (x,t) we have to check (|35|) at the time 

moment t — 0, i.e. to check the following equations: 

a k f a (x, 0) = a k f a (x, 0), a^f a (x, 0) = jjjgfe °- > , (37) 

Cat. 
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The first equation of (|37D follows from the definition of the function f a (x, t). To check the second equation 
of ( P?D we act on the function f a (x, 0) by operator 



ai : 



\m=0 



ml 



m=0 



ml 



and differentiate it with respect to variable a k 



df a (x,0) d 



E 



Oi 

m=0 \ i m=0 



da k \ — , , \/m\ 



E 



m k ip m (x,0) 



OL 



'ml 



L·* \/{m-e k )\ 



m=0 



where 



iV 



e fc = (0,... ,0,1,0,... ,0) 



(3Í 



(39) 



(40) 



(unit on the fc-th position). Comparing the two relations ( |38| ) and (^) we can conclude that the function 
f a (x,0) satisfies (|37|) . Making some transformations on both sides of (|3~7"D by means of the evolution 
operator 



Ü(t)a k Ü-\t)Ü(t)f a (x,0) = Ü(t)a k f a (x,0), 

df a {x,0) 



U(t)A+U-\t)U(t)f a (x,0) = U{t) 



dai 



(41) 
(42) 



and using the relation f a (x,t) = U(t)f a (x,0) we obtain ( pop . Therefore if the function f a (x,0) satisfies 
the function f a (x,t) satisfies (j35|) . Substituting the expression (|28|) for operator A(t) into the first 



equation of (^5|) we obtain 



À(t)f a (x,t) = -ihk p {t) ^^ +A x (t)xf(x,t) + S(t)f a (x,t) = af a (x,t). 



dx 



We can rewrite this equation in the form 

df a (x,t) 



-iU U ^ ' = A^ja - S(t) - A x (t)x)f a (x,t). 
It is easy to show that solution of equation (|44]) is given by the Gaussian function 

f a (x,t)=f 1 (a,t)expl -^A^^A^x + ^A^it^cL-Sit)) }. 



(43) 



(44) 



(45) 



To find the function fi(cc,t) we substitute this expression into the second equation of (p5|) 

a;(í) (-^^^) + K(t)*fi(<*,t) + s* (*)/!(€*,*) = + i(Aj(t))-Wi(«,í) 

Substituting the expression for —iU^- from equation ( f44|) we get the equation 



A*(í)A- i (í) a - <5(í) - A^í)* /i(a, í) + A^a^a, í) 



Sa 



(46) 



+ -(Aj(í))- i a; / 1 Kí), (47) 
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or 

= a;(í)V(*)(o - mh(<*,t) + (a:(í) - a;(í)V(*)a.(í) - ^(f)) -1 ) (48) 

Using the properties of matrices A x (t) and A p (í) from Appendix A one can show that the matrix in front 
of the x is equal to zero which can be inferred from the equality 

Aj»A:(í) - Aj^A^A^A^í) = Aj(í)A*(t) - A+(í)A x .(í) = l -E N . (49) 

Solving equation (|48|) we obtain the expression for the function /i(cc, í): 

(50) 



A(a,í) = exp j Va^A; 1 ^ + a r (ó*(t) - A^A; 1 ^)) + ¥>(*)} 



which contains the unknown function of time ip(t). Therefore, the wave function of multimode coherent 
state reads 

^«(aj,í) = exp<j -±-x T A^(t)A x (t)x + l -x T A p \t){ot - S(t)) + ^cx T A;(t)A^(t)a 

(51) 

+ a T (s*(t) - A* p (t)A;\t)ó(t)) - ±||a|| 2 + íp(t) 

To find the explicit dependence on time of the function (p(t) we substitute the expression (|5Ï"D into the 
Schròdinger equation ([!]). Using the following equalities (A and C are arbitrary symmetric matrices): 

pCp(xAx) = -2h 2 Tr(AC), 
p(xAx) = —2ihAx. 

we obtain the equation for the unknown function <p(t): 

^ ^(flW) - 5, p (í)) + iíÇ(í)V(t)í(í) - ^ < 5 T (í)(Aj(í))- 1 5 pp (í)A; 1 (í)5(t). 

(53) 



In order to solve this equation let us rewrite equations (15) in more detail, i.e. we get the evolution 
equations for the matrices 

À p (í) = A p (t)B xp (t) - A x (t)B„(t), À;(í) = A* p (t)B xp (t) - A* x (t)B pp (t), (54) 

À x (t) = A p (t)B xx (t) - A x (t)B px (t), A* x (t) = A* p (t)B xx (t) - A* x (t)B px (t), (55) 

and we get the evolution equations for the vectors 

6(t) = A p (t)c x (t) - A x (t)c p (t), 6*(t) = A* p (t)c x (t) - A* x (t)c p (t). (56) 



Using the evolution equation for the matrices (154]) we can transform the expression for which we calculate 
the trace of the matrix 

A-\t)A x (t)B pp (t) - B xp (t) = -A p \t)À p (t) (57) 
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Taking into account the known formula 

d 



Tr (A- 1 (í)i(í)) = - (ln det A(í)) , (5í 



we obtain the relation 

d 



Tr(A; 1 (í)A :i .(t) J B pp (í) - £^(í)) = --(lndet A p (t)). (59) 



To transform other terms on the right hand side of Ç|53f) we exclude the vector c p (t) and the matrix B pp (t) 



by means of (15) , i.e 



B„{t) = z^(aJ(í)A;(í) - A+(í)A p (í)), c p (í) = ili((í*(í)) r Af(í) - A+(f) (60) 

Substituting these expressions for the vector c p (í) and the matrix B pp {t) into the right hand side of equation 
(|53"D we obtain (thereafter we sometimes omit explicit dependence of some functions either on all or part 
of their arguments): 

+ -^(Ajr 1 (ajà; - a+à,) = * t a;a p - m - (í*) r « (ei) 

+ ^ t (a;a- - a;a-a A -), = - ^(^) + - (*T*) 



This expression gives the possibility to solve the equation fl53|) . Finally, we obtain the following expression 
for the normalized wave function of the coherent state 

if> a (x,t) = (27rn 2 )" Af/4 (det A p )- 1/2 exp j-^A;^ + ^/A^q - 5) + -a^A^a 

(62) 

+ c* T (> - A^Ç 1 *) - ~||a|| 2 + ^A^S - i||(5|| 2 + ijlm («jV) dr 

o 

We verify that the function ip a (x,t) is normalized in the partial case of the Hamiltonian without linear 
terms, i.e. we take the vector 6 = and calculate the integral 



(63) 



J \i/, a {x,t)\ 3 dx=\*{t)\*J expj-^^A.- (A* p ) 'A* x ] 

(A; 1 » - (A;)~V) + ^a T A* p A^a + ^(a*) T A p (A;) _1 a* - ||a|| 2 | dx 



where 



*(í) = (27rn 2 ) _Af/4 (detA p (í))- 1/2 . (64) 



To evaluate this integral, we transform the matrix of quadratic form in the exponent, using the properties 
of matrices A x and A p , i.e. 

A^A, - (A;)-^: = A^A. - A+tA+r 1 = A; 1 (A S A+ - A P A X + ) (A+)" 1 = -^-(A+Ap) -1 . (65) 
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Substituting (pq) into ( p3|) and evaluating the obtained integral we get 1 = 1. Thus the correct normal- 
ization of the wave function of the coherent state is proved. Using the formula for generating function for 
multivariable Hermite polynomials 



oxi> l —^a T Ra + a T Rx 



m=0 



ml 



(66) 



we can represent the coherent state (|62"D in the form (0). The coeficients in front of a m in (|32|) turn out 
to be the multivariable Hermite polynomials. Comparing the coefficients in (|66|) with coefficients in ([32]) 
we obtain the wave function of the Fock state 



ip n (x,t) 



x 



S-A p (A* p y l d*), 



(67) 



where ipo(x,t) is ip a (x,t) with ol = 0. 

We have shown that using time-dependent invariants method one can obtain the coherent and Fock 
states of the system with quadratic Hamiltonian. The problem of finding these states reduces to solving 
equations (EH) for the matrix íl(t) and the vector 6(t). As an example we consider the class of systems for 



which B px = B xp = and B pp and B xx are time- independent, permutable and non-singular. In this case 
for matrices A p and A x we get evolution equations with constant coefficients 



A, 



~A x Bpp -\- A.pB X p, 
- A. x Bp X -\- A.pB xxl 



or the system of equations in the form 



p 




B xp 




X 






- B px 



A T 

p 



Al 



Taking into account the initial conditions 

A p (0) = A p , 
the solution of this system of equations reads 



A*(0) = A a 



A?(*) 


= exp / 





Bpp 


'} 




A^(t) 




B xx 







Al 



(71) 



(72) 



(73) 



As shown in the Appendix B, the matrix exponent in right hand side of the equation (|73|) reads 

cos {^Bp P B xx t) - a/ BppB-^ sin ( yj B XX B PP t) 



exp 








Bpp 






J 












B xx 










B xx B p ^ sin ( \J Bp V B xx t) 



cos (y/B xx B pp t) 



(74) 



Thus one get the explicit form of the solution (|73| ) for the matrices 

A p (t) = Apcos [yjBp P B xx t^j - A x y/Bp P B- x 1 sin (y> ' B XX B VV , 
A x (t) = A p yJ B xx B~p sin (^J B PP B XX t j + A x cos Í^J B xx Bpp t 



(75) 
(76) 
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in the partial case of the quadratic multimode systems of the considered example. 

In this section we have shown that the wave function of the multimode coherent state of the 
quantum systems with quadratic Hamiltonian is given by Gaussian which depends on the parameters of 
symplectic transform which determine the linear invariants. We also have shown that the wave function 
(|6?D of Fick state of the multimode quadratic system is expressed in terms of multivariale Hermite poly- 
nomials which depend on the parameters of the symplectic transform, These two expressions for the wave 
functions of the coherent and Fock states are the main results of the section. 



3 Coherent and Fock States in the Probability Representation 

In this section we discuss the new probability or tomographic representation of quantum mechanics. We 
obtain the explicit expressions for tomograms of the coherent and Fock states of the quadratic mul- 
timode systems discussed in previous section in framework of Standard approach. In || an operator 
X = (Xi, Xpr), where n = 1, ...,N, is discussed for N = 1 as a genèric linear combination of position 
and momentum operators 



X-n /Ji n Xn ~\~ VnVr, 



(77) 



where \x n and v n are real parameters for all n, and X is Hermetian, hence observable. The physical 
meaning of the vectors jj, = (pi, ...,/ijv) and v = (v\, i>n) is that they describe an ensemble of rotated 
and scaled reference frames, in classical phase space, in which the position X may be measured. In |J it is 
shown that the quantum state of a system is completely determined if the classical probability distribution 
w(X , for the variable X is given in an ensemble of reference frames in the classical phase space 

(MDF). Such a function, also known as the marginal distribution function, belongs to a broad class of 
distributions which are determined as the Fourier transform of a characteristic function For the 

particular case of the variable (|77|) , considered in [|J, ||, 0, the scheme of [^] gives 



w(X,p, v) 



1 



(2ttH) 



N 



expi-ik 1 X) ( expiik 1 X) ) dk 



(7* 



where (A) = Tr(pA), p is the density operator. In [22| it was shown that, whenever X is observable, 
w(X, fi, v) is indeed a probability distribution. It is positive definite and satisfies the normalization 
condition 



J w(X, /x, u) dX = 1. 



(79) 



The defmition of the MDF allows us to express it in terms of the density matrix 



w(X,fj,, v) 



/ P(^>2')JJexp 

71=1 



flVr, 



dzdz'. (80) 



Recalling the relation between the Wigner function and the density matrix, the equation fl80D may be 
rewritten as a relation between w(X, /x, v) and the Wigner function, 



T exp[-ik n (X n - p n x n - u n p n ))W(x,p) 
n=i ( 27Th ) 



N 



51) 



Although the general class of distribution functions of the kind (|78|) was introduced, as a function of 



the density matrix, already by Cahill and Glauber in []22|, they did not analyze the possibility of a new 
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approach to quantum mechanics, in terms of such distribution functions, mainly because the invertibility 
of (^0) was not investigated. An important step in this direction is represented by [[J where Vogel and 
Risken have shown that for a particular choice of the parameters /i n and v n the marginal distribution 
completely determines the Wigner function. We make our previous statements more precise, saying that 
the MDF (the classical probability associated to the random variable X) contains the same information on 
a quantum system as the density matrix. Marginal distribution w (X, /x, v) can be also expressed through 
wave function of the system. This expression is given by 



w{X,n, u,t) 



1 



(27rh) N \u 1 ■ . . . ■ u N \ 



i/>{y, t) exp<| l^y T MN' l y - ^N^X \ dy 



52) 



where the diagonal matrices which contain the parameters determining the reference frames in the phase 
space of the multimode quadratic system 



M = diag{/ii, . . . ,/íjv}, 
N = diag-{>i, . . . , v N }. 



i3) 
54) 



are introduced. Looking at the genèric tomogram (|8~2|) it is obvious that to evaluate the tomograms 
corresponding to coherent state ^ a (x, t) and the Fock state ip n {x, t) we need to evaluate the following two 
integrals: 



and 



I a (t) 



I n (t) 



> < a (y.t)ow{ ^y T MN- l y - \y T N~ l x \ dy 



i> n (y, t) exp<| ^y T MN- l y - ^fN^X \ dy. 



15) 



16) 



Remembering the expression for the generating function ( p2[) we have the relation 

00 n 



Tl=0 



57) 



which shows that the integral (|85| ) is the generating function for the integrals (|36|) . In view of 
obtain expression for in terms of derivatives of I a 



we can 



d 



where 



d 



dct 7 



d 



a=0 



(e^ M / Q (í)). 



d 



da n N » 



dcx n da^ 1 

According to ( |62|) the wave function of multimode coherent state ip a (x,t) can be represented as 

V»«(aJ,í) = y(t)exp^-^x T A x + ^x T A p \a - S)+K(a,t)\, 
where we introduced the matrix 



19) 



(90) 



A = iA p 1 A X 



(91) 
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and the function 



K(a,t) = l -OL T h* p k p l cx + a T (S* - A^S) + ^6 T A* p A; l 6 

t 

1 "~ '|2 , • f r f^sA j 1 u ||2 



, )ll S\\ + i Im d* ) dr - -\\a\\. (92) 



which does not contain the dependence on JT-variable. Using the known formula for multidimensional 
Gaussian integral 



exp>{ — -a: Ae + x a > dx = y ^— exp< -a A a}>, (93) 



we can evaluate the integral fl8~5]) . After some transformations using properties of the matrices A p and A x 
we can obtain the following result: 



where we introduce the function 

E(X, fi, v, a, t) = exp j-^<5(X, [i, v, et, t) + K{ot, t) \, (95) 
which is determined by the function 

Q(X,fj,,u,a,t) = -\{X- NA p 1 (a-6)) T N- 1 (A - ï MN- l y l N^ 1 (X — NA p 1 (a. — 6)) . (96) 
For the tomogram of the multimode coherent state w a (X, /x, v,t), we get the expression 



w a (X,fj,,u,t) = ^ n , ^J^^' \E{X,fj.,is,a,t)\ 2 , (97) 

(2-iïh 2 ) ' \v x ■ ■ ■ ■ ■ u N \ -y/detc+(/Lt, i/,í)n(/x, u,t) 



where the matrix H(/x, u, t) is given by 

E{fi,u,t) = iA x (t)N -iA p {t)M. (98) 
On the other hand, we can represent w a (X, //, u, t) as multidimensional Gaussian distribution function, 



l.C 

1 



(27r) JV detS 



w a (X, fi, u, t) = - = exp<j -1(X - X ) T S- 1 (X -X )}, (99) 



2 



where the dispertion matrix depends on the reference frame parameters and the parameters of symplectic 
transform 

£0, u, t) = £ 2 ~ + (/x, v, t)E(fi, u, t), (100) 
and the mean value of the multidimensional random variable 

X (ii, v, ol, t) = HZ+il·L, v, í)S(/u, u, t) (S- 1 ^, v, t)(a - 6(t)) + (S*)" 1 ^, i/, t)(a* - ó*(t))) . (101) 
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is determined also be the symplectic transform parameters and depends on the reference frame parame- 
ters. Using expression for the tomogram w a (X , /J,,u,t) and formula for the integral fl88|) we can obtain 
tomograms of the Fock state: 



w n (X, n, v, t) = w (X, u, v, t)—r 



ni 



h 



(102) 



where Wq(X, u, u, t) is w a (X, u, u, t) with a = 0. For one-dimensional case this formula takes the form 

2 



w n (X, n, v, t) = w (X, /i, v, t) 



1 



2 n n! 



n \hy/2\Ç\ y/2\t\ 



(103) 



where the matrix H, is the one-dimensional matrix, and it is denoted as £. The numbers £ and S depend 
on time. 

The explicit expression of the tomogram of the Fock state of the multimode quadratic systems in terms 
of the multivariable Hermite polynomials ( |102| ) is the main result of this section. 



4 New relations for Hermite Polynomials 



In this Section we get some new formulas for multivariable Hermite polynomials using results obtained in 
previous Sections. In section 1, considering the problem of finding coherent and Fock states, we have found 
these states for matrices A p and A x which satisfy conditions (0) and ( 168 ). If we chose other matrices 
A' p and A' x which also satisfy the same conditions, we obtain other coherent states and other Fock states. 
We want to find out how "new" Fock state ip' n with matrices A' p and A' x can be expressed through "old" 
Fock states tpm with matrices A p and A x . We consider the case c = (hence 6 = 0). 

Making symplectic transformation of annihilation and creation operators a and a + of the form 




S 



a 



S 



S p 




Q* 


Q* 



(104) 



we obtain "new" operators a' and a' . These operators must satisfy condition (|Ï8[) and therefore the 
matrix S must satisfy the following condition: 



SY>2nS T 



(105) 



or 



S P S X 



c çT Q* c+ c* c+ 

>Jx>Jp — àp^x D x°p 



Ç Q+ C C + 
OpOp °x>->x 



n* qT C* C J 

OpOp o x o x 



E 



N- 



(106) 
(107) 



It is easy to show that iV x iV-matrix S p is non-singular. In fact, let z ^ is the vector for which SpZ = 0. 
It means that the matrix S p has no inverse, i.e. it is singular one. In this case equation ( |107|) gives for 
scalar product of vectors the inequality 



(2, SpSpZ^j (z, S x S x z) — [S£ z, S x z) — \\S x z\ 



> 0. 



108) 



This inequality contradicts to the inequality — 11*9+2 1| 2 ^ 0. Thus, we proved the nonsingularity of the 



matrix S p . 
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The functions ip n (x, t) form complete set of function and we can represent the function ip' n {x, t) as the 
series 



il>' n (x,t) = c nm i! m (x,t). 



(109) 



m=0 



The complex coefficients c nm have the physical meaning of transition amplitude from an initial Fock state 
\n)' to a final Fock state \m). According to (|67|) and ( |104j ) "new" Fock state ij)' n {x,t) is given explicitly 
by the formula 

i(/ n (x,t) = (27rft 2 ) _JV/4 (det (S p A p + S x A* p ))~ 1/2 exp \-±_ x T (S p A p + S^;)' 1 (S p A x + S X A* X ) x 



1_ {-(s;a* p +s* x a p )(s p a p +s x a* p ) f i 



(ApSp + ApS x ) x j . 



h 



;no) 



As it is known the decomposition coefficients c nm are given by scalar product of wave functions if) m (x,t) 
and ip' n (x, t) which is expressed in terms of the overlap integral of two multivariable Hermite polynomials 
and Gaussian of the form 



c nm = (^ m \^' n ) = {2nh 2 ) N/2 (det(S p A p + S x A;))' 1/2 (detA;) 1 



1/2, 



1/2 1 1 



nl vml 



exp j —^x T (-(A* p )~ 1 A* x + (S P A P + SzA;)- 1 (S p A x + S x A x ))x 



{-(s;a;+s*a p )(s p a p +s x a;) l } í i 



x H, 



Here we used Dirac notations 



h 



{(f | ip) = / (p*(x)ij)(x) dx 



(111) 



(112) 



for the Standard scalar product of two vectors </?), |^) G L2(C N ) written in position representation. 
Evaluating this integral we obtain for the transition amplitude 



yrúrrú 



[t n s m ] exp<j ^t T s* x s;H + a T s;H - \s T S; l S x s 



^=7= ^(0,0), (113) 



yj det S'p 
where the 2iV x 27V-matrix 

0*0-1 ( qt\~ 1 

~ J x J p \ J p ) 

Q—l Q—X C 

°p °p °x 

is expressed in terms of four block-matrices. The normalization condition of ip' n (x,t) is reduced to the 



(114) 



equality Yl \c nm \ 2 
the following from: 



1. Using the expression (|113|) for coefficients c nm we can represent this equality in 



j +00 
|detS p | = - 



<L)(0,0) 



m=0 



ml 



(115) 
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This form gives a sum rule for multivariable Hermite polynomials. Let us consider the partial case n = 0. 
From ( |113| ) we obtain transition amplitude from ground state |0) to an excited state \m). 



COr 



*/ det S p 



ml 



(116) 



The expansion 



;ii7) 



takes the explícit form 



í 1 



m=0 



-1 C C /iT\ — 1 



expj— x T {S p A p + S^)" 1 ^^) * ) = yjdet (E N + S x A* p A^S^) 



2h 2 



'118) 



m>0 



ml 



x 



The obtained result provides the partial case of sum rule for multivariable Hermite polynomials. Normal- 
ization condition for the new ground state of the quadratic system ip' (x,t) is equivalent to the equality 



|det S p \ = 



hL Sp ~ 1Sx \o) 



m=Q 



ml 



(119) 



which is new sum rule for the Hermite polynomials. We can verify that this equality is true in some 
particular case. For example, if N = 1 and the matrix 



S 



cosh 9 sinh 9 
sinh 9 cosh 9 



(120) 



i.e., S p = cosh#, S x = sinhó 1 , the Hermite polynomials reads 



H {^o }{0) = ^^y /2 Hi 



(0) 



\m (2m)! 
' 2 m m! 



n = 2m + 1 
(tanh#) m n = 2m 



(121) 



and we obtain from ( |119|) the decomposition of the form 



coshfl = 2 ,l 2m) l 2 (tanhg) 2m 



m=0 



2 2m (m\y 



(122) 



Taking into account that cosh 9 — 1/ y 1 — (tanh 9) we can conclude that ( |122| ) is true, since expansion 
of function l/\/l — x is given by the obvious relation 

I ^ (2m)\ 



+oo 

E 



y/l-x ^-1 2 2m (m!) 

m=0 V / 



(123) 
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Let us consider the case n ^ 0. Using the relation between two-dimensional Hermite polynomials and 
Legendre polynomials: 



^ } (0,0) = (min(n,m))!(-l)^/Vr/ 2 r 2 m 2 /2 f ~4 



\ r ll r 22 



(n+m)/4 



P, 



\n-m\/2 



7*12 



(n+m)/2 ^ r 2 2 _ rilí . 22y / 



1 , (124) 



where the 2 x 2-matrix R has the form 



i? = 




ris 




r 2 \ 


r 22 



we obtain for matrix S given by (|120 ) the expression 

<L)(0,0) = (min(n,m))!(-l) 5 «/ 4+3 -/ 2 P ( l: + -:J/ 2 

where the matrix elements in (|125|) are taken in the form 

— tanh 9 



(125) 



(126) 



R 



i 



cosh 9 



cosh d 

tanh# 



(127) 



Comparing the expression for transition probabilities obtained in probability representat ion of quantum 
mechanics and in usual representation we can find other new relation for Hermite polynomials. Transition 
probability in usual representation is given by the known relation 



w nm = \ (^ n (ti) | ^m(*ü))| 



128) 



Therefore to represent the transition probability in explicit form we need to calculate the following overlap 
integral of two wave functions: 



(V'n.(íi) | ïpmih)) = / ip n {x,t 1 )ip rn (x,t 2 ) dx 



z/í) detD 

+ l-S* T (t 2 ) E*D~ 1 6* (ix) - U T (tx) D- X S* (t 2 ) 



= exp(V (tjD^ES (h) 



(129) 



|<5(í 2 )|| 2 + i J Im(ïV) dr\ 



HÏ R} Ju,v) 



t-2 



// j T (n,m) V ' 
yn\m\ 



where we introduced new 2N x 2iV-matrix 



R = R(ti,t 2 ) 



D- l E 

Ud t )~ 1 ew- 1 



(130) 



and the 2iV-vector 



v) \v(ti,t 2 ) 



-E+iD+y 1 -Ud*)- 1 

1 1 



-E{D*y l 



(d-^EÒ (h)-6* ( tl )-~D- l d* (t 2 ) 



E*D- 1 6* (t 2 )-S(t 2 )--D- 1 S(t 1 ) 

n 



(131) 
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The N x A^-matrices 



d = D{t u t 2 ) = a; (í 2 ) Ki {t x ) - a; (í 2 ) (í x ) 



(132) 



£ = JSfa, í 2 ) = A; (í 2 ) A+ (iO - K* x (í 2 ) AJ (Íí) . (133) 

Are used to define the 2 N- vector (|131|) , and these matrices are connected with symplectic matrix deter- 
mining the linear integrals of motion of the system under consideration. On the other hand, the transition 
probability from the initial state n) to the final state |m) in probability representation is given by the 



probability from the initial state n) 
expression 

1 J w n (X', n, u, ix) w m (X", -n, -u, t 2 ) e i{x ' +xll) dfi du dX' dX 



(2tt/ 

where we used the notation 



(134) 



TV 

n 

k=l 



e Xk = e k = x 



N 

E x k 



(135) 



Using expression ( |102|) obtained for the tomogram of the Fock state w n (X, fi, u) and substituting it into 
( |134[) , we obtain the equality 

/|iíi (HT(íl)rlH+(<l)} Q(s + (t 1 ))"^ / +(s + (í 1 ))- 1 (s + (t 1 )^(íi) + ^(í 1 )^(t 1 

i (ST(Í2)rV(Í2)} f-i(H+( Í2 ))- 1 X"+(H+( Í2 ))- 1 (5+(í 2 )<5( Í2 )+^( Í2 )r(í 2 

V h v ' 7 (136) 



x Wo (X', /x, i/, íi) Wo (X", -M, í 2 ) e- l(x ' +x " ) du du dX' dX" = ) 

n iy det D 



x 



For the partial case 6 = 0, this general formula gives 



-{(^)) 



X W 



(X', /x, í/, ti) w (X", -jí, -i/, £ 2 ) e í(x ' +x " ) rf/x di/ dX' dX" 

1 1 2 



(137) 



<L)(0,0) 



The obtained two formulas are new mathematical results which are obvious from the point of view of 
quantum transition consideration. But purely mathematical derivation of the obtained formulas for mul- 
tivariable Hermite polynomials is not that easy problem. 

5 Examples of Systems with Quadratic Hamiltonian 
5.1 One-dimensional Harmònic Oscillator 



In this section, we consider some particular cases of systems with quadratic Hamiltonian, such as harmònic 
oscillator, particle in elèctric field. We calculate tomograms of coherent and Fock states of the system. 
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Our first example is a harmònic oscillator with time-dependent frequency and driving force. This 
oscillator is described by the nonstationary Hamiltonian 

5(0 = ^ + =^ + /(«)£■ (138) 

This Hamiltonian is used to describe trapped ion. For the Hamiltonian the matrices A p used in previous 
sections and A x are given by the expressions 

A p = -^==e, (139) 
V zmuúh 

A x = —j=L=mè, (140) 

where the function of time e satisfies the equation 

i + cü 2 e = 0. (141) 

The initial conditions are taken in the form 

e (0) = 1, è(0) = cu(0)i. (142) 
For the one-dimensional matrix H, dispersion matrix £ and vector Xq, we have 

'mèl·' + efi), (143) 



h 

2mw 



S = 7T^— |mèí/ + e/i| 2 , (144) 



X = \/— ? f(méV + eV)(«-í) + (rraèz/ + efi) (a* -5*)). (145) 

V 2muj \mèv + e/il V ' 



In this particular case all the matrices and vectors are one-dimensional objects. The wave functions of 
coherent states and Fock states functions are given in explicit form 



.Jmuj 1 (i mé 2 2muj x(a — 5) 1 e* 2 



t 

+ a(s* + Çs^J - ~5 2 - ^\5\ 2 + ijlm (<W*) dr - i|a| 2 J, (146) 

o 



. miü lli e* \ I i mé 2 l2mwxó 
Vn{x,t) = J ——=—==\ \ — \ exp<^ — — x 



irh y^e y / 2 n n\ \ V £/ (2fi £ \ h e 



t 



1 * > - + i í Im fe*) drW (M* + 5 Jl±ll ) . (147) 



2 e 2' 1 i V ) ) n \\ h \e\ y/2\ 
o 

As in general case the coherent state is described by Gaussian wave function. The Fock state is described 
by one-dimensional Hermite polynomial. The parameters of the wave functions depend on the property 
of the linear integrals of motion of the oscillator. Substituting the expressions ( |143| )-( |I~4"5| ) into fl9"9|) and 
(|102[) , we obtain the tomogram of the coherent state 



w a (X,n,u,t) = J^ . .\ , exp \~— AX-X T} (] IS) 

nn \mev + efi\ { fi \meu + efi\ 



18 



and the tomogram of the Fock state 

1 



w n (X, fi, v, t) = w (X, fi, v, t) 



2 n n\ 



mu X (mè*i> + e*fi)8 + {mèv + efi)8* 



h \mèis + efi\ 



\mèis + ejx\ 



(149) 



of parametric oscillator under consideration. 

If driving force is absent then 5 = and we have the following expressions instead of previous ones. 
The wave function of the coherent state is given by the Gaussian 



mu 1 
nh -Je 



i mè 2 / 2mw xa 1 e* 2 1 . 



a \a\ 

he 2 e 2 1 1 



(150) 



The wave function of the Fock state is expressed in terms of Hermite polynomial 



moo 1 1 



cxp 



% me 

2h~7 



x 2 \H n 



muj x 
h \e\ 



nh ^fè V2 n n! V V e 

The tomogram of the coherent state coincides with Gaussian distribution function of the form 

1 



(151) 



w a (X,fi,u,t) = 



mu 



■ cxp 



mu 



nh \mèv + e/i 

The tomogram of the Fock state of the parametric oscillator reads 

1 



h \mèu + eji\ 



(152) 



w n (X, n, v, t) = w (X, h, is, t) 



2 n n\ 



H„ 



mu X 
h \mèv + e/i\ 



(153) 



For the case of constant frequency u = const ^ 0, we have e = e tuJt and the expressions for the wave 
functions of the coherent and Fock states and the tomograms of these states are given by simple formulas. 
The wave function of the coherent state of the harmònic oscillator in position representation reads 



mu 



ip a (x,t) = (/— -exp<^ — X + 



nh 



mu 



2h 



2mu 



e~ iuJt xa - -e-^a* - -\a\' - -iut >. 



-2iort ,2 



h 



(154) 



The wave function of the stationary state of the harmònic oscillator describing for the one-mode electro- 
magnètic field the n-photon state in the position representation has the Standard form 



Ís n {x,t) 



mu 1 



-iüjí(n+i) 



exp 



mu 
'^2h 



x 2 \H„ 



nh y/2 n 7l\ 

The tomogram of the coherent state of the harmònic oscillator reads 

1 



mu 



-x 



w a (X,n, u,t) = 



mu 



nh 



' ji z + m?u 2 is 2 



cxp 



mu 



h /i 2 + m?u 1 v 2 



h 



(X-X ) 



(155) 



(156) 



The tomogram of the n-th excited state of the harmònic oscillator takes the form 



w n (X, fi, is, t) = w (X, fi, is, t) 



2 n n\ 



H„ 



mu 



X 



h 'fi 2 + m 2 u 2 v 2 



(157) 



Thus we obtained explicit forms of tomograms for the specific quantum states of parametric and harmònic 
oscillators. 
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5.2 Charged Particle in Nonstationary Electric Field 



Important partial case of quadratic system is the problem of motion of a charged particle in a homogeneous 
time-dependent elèctric field. The Hamiltionian of the particle reads 

^2 

(158) 



P 



H(t) = ^ + F{t)x, 
2m 

where F(t) is nonstationary elèctric field. Here the mtarices and vectors determining the genèric Hamil- 
tonian of a quadratic system (||) are reduced to the numbers B pp — 1/m, B px = B xp = B xx = 0, c p = 0, 
c x (t) = F(t) and for matrices A p , A x and vector 5 we have 



A x {t) = A x , 



A P (í) = A x t + A p , 

m 



(159) 
(160) 



6{t) 



A n -—t) F(t)<It. 



m 



(161) 



Using expressions ( |159| ) one can obtain the wave function of the coherent state of the charged particle 



ip a (x,t) 



1 



+ oc | 5* - 



exp 



A p -^t 

P m 



A* - 

_P m S 

A„-^t 

P m 



A, 



2hA v - à*t 



lA*-^t „ 

P m x2 



£ + - 



i x(a — 5) 



+ - 



P m ^ 2 



hJL-M 2A v -^t 

P m v m 



OC 



2A p -^t 



5 2 --\5\ z + i / Im 55* dr--\a 



(162) 



and the wave function of the Fock state of the particle 
1 1 1 



ip n (x,t) 



A* - &t 

V m 



V P m 



A p -^t 

P m 



P m 







x H„ 



x 



hy/2\A„- ^t\ 

v \ P m 



+ i 



A* A % 
P m 



exp 



A, 



x5 



2hA v -^t 

P m 



X — — - 



P in 



V2\A p -%t\ 



(163) 



These states are nonstationary normalized states of the charge in the elèctric field. Tomograms of coherent 
and Fock states are given by the formulas fl9~9|) and ( |103| ) in which the matrix S, the dispersion matrix E 
and vector Xq take the following vàlues: 



iA v v — i ( A p — 



A x 



tí 



A x v 



A, 



m 

Ax 
m 



t //, 



t /i 



X = h ((iA x v - i Ía p - ^-t\ (a-5)+ (iA x v -i(^A p - àltj fij 



oc 



(164) 
(165) 
(166) 



The parameters A x and A p determine the degree of squeezing in the coherent state of the charge in the 
elèctric field. 
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6 Conclusion 



We have shown that for multimode systems with Hamiltonian which is quadratic in position and momentum 
operators the quantum states can be described by symplectic tomograms. For coherent states these 
tomograms are the multivariable Gaussian distribution functions. For the Fock states the tomograms 
are expressed in terms of multivariable Hermite polynomials. New results of our work are the explicit 
formulas for tomograms of coherent and Fock states of multimode nonstationary quadratic systems and 
new relations for multivariable Hermite polynomials. Since many physical systems including charge moving 
in elèctric and magnètic fields and photons in a resonator with moving boundaries are described by 
quadratic Hamiltonians the results of this work can be applied to consider the evolution of quantum states 
of these physical systems in the tomographic probability representation. The mathematical formalism 
of these quadratic systems is based on properties of multivariable Hermite polynomials (p3|, p^j , p5| . 



2(|). Due to this the new formulas for multivariable Hermite polynomials obtained in the tomographic 
representation can be used in analysis of the different physical systems with quadratic Hamiltonians. It 
is worthy to consider some other particular cases of multimode quadratic systems. We hope to do this in 
forthcoming paper. 

A Properties of symplectic A-matrices 

Using definition of operators a and a + we obtain the following properties of matrices A p and A x : 

A x A T p - A P A T X = 0, (167) 



A X A+ - A p At = ~E N . (168) 



Using these two properties, one can prové that the matrices A p and A x determine the real symplectic 
transform of position and momentum, i.e. 

(i) Matrices A p and A x are non-singular, 

(ii) ApA* = A+Ap, A T X A* X = A+A x , 

(iii) A X A P — A X A* = A p A x — A p A x = jEn, 

We prové first the property (i). Suppose that matrix A p is singular. Then the transposed matrix A^ 
is also singular and there exists a nonzero vector z such that A^z = 0. From ( 168|) we have the chain of 
the relations for the scalar products 

+ ^\\z\\ 2 = ~ (z,z) = {z,A x A; Z ) - (z, A p A+z) = - {A;z,A+z) = (169) 

These contradictory relations show the nonsingularity of the matrix A p . The proof that the matrix A x is 
also nonsingular is analogous. 

Let us prové the property (ii). In view of nonsingularity of the matrices A p and A x we can write ( |168j ) 
in form 

A-t ~ A p 1 A x A p l· = -A p 1 . (170) 



From f [L67|) it follows that 



A p 1 A x — A x Ap . (171) 
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Also we can obtain that 



A + _ 4*-l A* A + 
A x — A p A x A p ■ 



Substituting this expression for and the equality (|171|) in ( |170| ) we get 
This formula can be rewritten in the form 

A*-l I A* AT A*AT\ AT' 1 A+ _ i A-l 

A v \ A x A P A p A x ) A P A p — j^p ■ 

Using hermitian conjugate of (|168|) we see that the above equation gives the relation 

A*^ 1 4 T_1 A + = A^ 1 
p p p p 

or 

4 + 4 __ at a* 

The second equality in (ii) can be proved analogously. 

Let us prové the property (iii). Multipling ( |168| ) by from the left side we get 

AlçA x Ap — A+A P A+ = —jA^. 

Using the property (ii) we get the relation 



or using ( |167| ) we get 



4 + 4 4+ _ A T A* A + - -A + 

A x ApA x A x A x A p — ^xi 



A + A A + — A T A* A + — — A + 
j-í x -n-p-n-x ^-x^-p^-x fi 



Since the matrix A x and consequently the matrix A+ are nonsingular we obtain the relation (iii). 
Calculating the time derivative of the product of the matrices A(t)E 2 ArA T (í) we get the equality 

j t (A(Í)£ 2A A T (Í)) = A(í)£ 2iV A r (í) + A(í)£ 27V À r (í) = 0. 

Taking into account the initial condition for matrix A(í) ([ÜJ), we obtain 

A(Í)£ 27V A T (Í) = E 2N . 

Making the same procedure with the product of matrices í2(t)S 2 jyí2 T (í) we obtain 

i 



Q(t)E 2N Q T (t) = -E 2N . 



Rewriting (|182|) in more detail, we have 



A X A^ - A P A^ = 0, 



A^AJ - A p A+ = -~E N . 



172) 

173) 

174) 

175) 
176) 

177) 
178) 
179) 



180) 
181) 
182) 



183) 
184) 



K 

These equations are similar to ( p.67| ) and ( |168| ) and therefore the matrices A p and A x possess the same 
prop ert ies as A p and A x , i.e., 

(i) Matrices A p and A x are non-singular, 

(ii) aJa; = a;a p , A t x A* x = A+A x , 

(ííí) AtAp - aIa; = A T p Ai - a;a x = 1 -e n . 
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B Calculation of the Matrix Exponent 

In this appendix we calculate the exponent 





• 


—B pp 




+00 1 

n=0 


—Bpp 


n 


Epp{t) E px (t) 


exp < 






■! 




f 1 = 








B xx 




B xx 




E xp (t) E xx {t) 



To calculate this exponent it is worthy to note that 



—B 



B xx 

—B. 



2k 

PP 
2k+l 



~(BppB x 



( B xx Bpp J 



pp 



B xx 



[—Bpp) (—B XX B, 



pp) 



~B XX ) ( BppB x 



In view of this we get 
1 
2! 



]_ 1 2 ( 

E pp (t) = E N — —Bp P B xx t 2 + —(Bp P B xx ) t 4 — . . . = cos yy/ B PP B XX t 

E px (t) = —Bppt + -^BppB xx B pp t 3 — -^B P p(B xx Bp p yt 5 + . . . = — yj B PP B~^ sin ^ y/ B XX B, 



E X p{t) B xx t ^B xx BppB xx t -\- ^B xx (BppB xx ) t -{-... 



B xx Bp V x sin ( yj B VV B XX t 



E xx (t) = E N — —B xx B P pt 2 + —{B xx Bppft 4 — . . . = cos {^J B xx B pp tj 
Thus for the matrix exponent (|185|) we have the expression 



exp 

Solution (|73|) is 



- 


Bpp 


B xx 





is 













cos (y/B pp B xx t) - y/ BppBj sin ( y/ B XX B PP t) 

y/B xx B-p X sin (y/Bp P B xx t) cos (y/ B xx B pp t) 



or 



cos (y/B pp B xx t) - y/ BppB-i sin (y/B xx B pp t) 

y/B^B-p 1 sin (y/Bp P B xx t) cos (y/B xx B pp t) 

A p (t) = Apcos [y/Bp P B xx tj - A x y/BppB~l sin [y/ B XX B PP t 
A x (t) = Ap^B xx B-p X sin (^y/B pp B xx tj + A x cos (a/ B xx Bpp t 





Al 




Al 



If we take A p and A x as 



we have 



A,, 



i / ^pp ^ra j 



.4, 



y[~Bp~pB 



1 

11 > 



A p (t) 



. 27i 

These expressions were used in the main text. 



-^={/Bp P Bj exp |za/B^B^íJ , 
A^í) = -j=^B~ p l B xx exp {i^/B pp B xx t} . 
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